Dynamics after a sweep through a quantum critical point 
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The coherent quantum evolution of a one-dimensional many-particle system after slowly sweeping the Hamil- 
tonian through a critical point is studied using a generalized quantum Ising model containing both integrable 
and non-integrable regimes. It is known from previous work that universal power laws of the sweep rate ap- 
pear in such quantities as the mean number of excitations created by the sweep. Several other phenomena are 
found that are not reflected by such averages: there are two scaling regimes of the entanglement entropy and a 
relaxation that is power-law in time rather than exponential. The final state of evolution after the quench is not 
characterized by any effective temperature, and the Loschmidt echo converges algebraically for long times, with 
cusplike singularities in the integrable case that are dynamically broadened by nonintegrable perturbations. 
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A many-particle quantum system evolving at zero temper- 
ature can demonstrate various forms of approach to equilib- 
rium even with no loss of phase coherence. This phenomenon 
has been studied in most detail experimentally and the- 
oretically El H S H 01 f° r systems prepared by a quan- 
tum quench across a phase transition. A system is prepared in 
the ground state for certain parameter values, which are then 
rapidly changed to values for which the ground state is in a 
different phase. Ultracold atomic systems are especially valu- 
able for experiments in this area because they can be treated as 
closed quantum systems on rather long time scales compared 
to the basic dynamical time scales of the system. Two basic 
notions in the literature are that many properties equilibrate af- 
ter the quench exponentially in time and that the system ther- 
malizes (the final state can be described by an effective tem- 
perature). The notion of thermalization is generally not a well 
defined one for integrable systems and only non-integrable 
systems in the thermodynamic limit can be described by an 
effective temperature after a quench @|. Non-integrability 
alone may not be a sufficient condition for thermalization in 
finite systems 

This paper concerns the coherent zero-temperature dynam- 
ics of states prepared through a different process that has also 
attracted recent attention EH EH 0, Q: rather than 
quenching instantly across a quantum phase transition, the 
Hamiltonian is changed smoothly across the transition at a 
constant rate V. For a second-order transition (critical point), 
various averaged physical quantities such as the excitation 
density and energy show power laws in rate T as V — > 0, with 
exponents determined by the quantum critical point's univer- 
sal physics. The system evolves after such a sweep to a steady 
state for some quantities, such as the Loschmidt echo defined 
below, but its energy distribution remains non-thermal. The 
main result of this Letter, using a generalized quantum Ising 
model as an example, is that many features of the resulting 
evolution including the approach to equilibrium differ from 
the case of an instant quench and are not captured by the 
simple averaged quantities studied previously. Integrable and 
non-integrable systems evolve differently and our results sug- 
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FIG. 1: The panel shows the energy levels e n (g) and gaps of eigen- 
states \ipn(g)). The energy density Eb(g,T) of a wave function 
\ip'o(g,T)) resulting from a linear sweep at a constant rate T and 
the ground state energy density Eo(g) are shown in red. 



gest a sharp dynamical probe of these differences. 

We consider a system with a tunable parameter g that be- 
comes critical at a value g c and study the dynamics as the 
Hamiltonian is swept from one side of the critical point to the 
other. A particular example of such a system is the generalized 
quantum Ising model in a field described by the Hamiltonian 



H = J 



N 

£ 



g (cos (pa i + sin (pa\ 



(1) 



The system is critical at g = g c = for any value of <f>. For 
= 0, the model Hamiltonian is the Ising model in a purely 
transverse field (the "quantum Ising model") which is inte- 
grable. For any 0^0, the lattice model is no longer inte- 
grable (except at the critical point where g = 0). We fix the 
value of <f> and do not write the <j> dependence of the quantities 
explicitly in the discussion below. The correlation length goes 
as £(ff) ~ \g-9c\~ v and the gap as A(g) ~ ~ \g-g c \ zv , 
close to the critical point, where z is the dynamical exponent 
and v the correlation length exponent (there are two values of 
v depending on <f>, v = 1 for <f> = 0, ir and v = 8/15 other- 
wise). At each g, the system has a ground state |-0o(<?)) with 
energy eo(g) and excited states IVVi (<?)); the lowest excited 
state \ipi(g)) has energy A(g) (FIG.[T|i- 

We study the dynamics of energy, entanglement entropy, 
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FIG. 2: The difference AE between the energy density E' (gf,T) 
of the final state after sweeping with a rate F through the critical 
point and the ground state energy density £o (<?/). The dashed lines 
show the expected asymptotic behaviors AE ~ VT (<f> = 0) and 
AE ~ p 8 / 23 (0 > resulting from different correlation length 
exponents in the "thermal" and "magnetic" directions. 



and wavefunction overlap during and after an adiabatic sweep 
across a quantum critical point. We derive several analytic 
expressions and compare to numerical results for a sweep 
through the critical point of the Hamiltonian (Q]). We assume 
that the parameter g depends on time t as g t = g% — Tt/J, 
where <?i is some initial value of g and % = 1, until some 
final value g f is reached, after which the Hamiltonian is con- 
stant for some "wait period" . We are interested in the adia- 
batic limit r — > 0. For the numerical calculations, we use the 
recently introduced infinite Time-Evolving Block Decimation 
(iTEBD) algorithm 11511 . This method uses variational wave- 
functions based on matrix product states (MPS) and exploits 
translational invariance for efficient simulation of infinite ID 
systems. Its errors result from finite entanglement rather than 
finite size |3 17]. Away from a critical point, the entan- 
glement entropy of the exact wave function is finite, and the 
iTEBD algorithm becomes very accurate. 

The expected energy density E' (g tl T) of \ip' (g t , T)) at a 
time t (i.e., the expectation value of the Hamiltonian with 
coupling g t in the time-evolved state) when the system has 
crossed the critical point is a natural quantity to consider. We 
use a sca ling relation for slow sweeps in this model derived 
in Refs. Il2l 1 1311 . Let the system be in the state \ipo {gi)) at 
t = and change g t linearly with rate V — > 0. Whilst g t is 
on the same side of the critical point as gi, the system will 
be in state \ipo(gt)) with exponential accuracy as guaranteed 
by the adiabatic theorem. However, adiabaticity breaks down 
at the critical point since the gap vanishes, and once g t is on 
the other side, the system is in a state \i[;' (gt, T)). Thus, ex- 
citations (well-defined in the = model) are created in the 
system upon crossing the critical point. The number of ex- 
citations is ?i ex = CT du ^ zl/+1 '> lUl], where d is the spatial 
dimension. Note that the scaling parameter in the above equa- 
tion is the rate T and not \g — g c \ since n ex is non-zero only as 
a consequence of sweeping across the critical point at a finite 
rate. Since the excitations in the transverse Ising model can be 
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FIG. 3: (a) The entanglement entropy as a function of the inverse 
rate. The dashed line in the upper panel shows the expected asymp- 
totic behavior S = — 1/24 ln(r) + const, (b) Entanglement entropy 
as a function of time. The dashed lines indicate the time at which the 
final value of (gj — —0.4) has been reached, and the Hamiltonian 
remains unchanged thereafter. 



interpreted as free domain wall excitations, the energy density 
of the final state is expected to be proportional to n ex , so 

AE = E' {g t ,T) - E (g t ) ~ n ex A( 5t ) ~ T d ^ v +^A(g t ). 

(2) 

The gap A(g t ) depends only on the instantaneous value of gt 
and does not scale with T. Eqn. (O follows from noting that as 
T — > 0, there are a small number of excitations that exist only 
in a small band of states of vanishing width above (<?))■ 
Numerical results are shown in FIG. |2] The energy difference 
AE between the actual ground state and the state after the 
sweep for the transverse Ising model (<f) = 0) is in a good 
approximation proportional to vT. The critical exponents for 
the transverse Ising model are z = 1 and v = 1 and n ex oc 
VT (see also Refs. H 1 lL Il2ll ). Thus the numerical results are 
consistent with the scaling in Eqn. (0. 

We now consider the entanglement entropy between the left 
and right halves of the infinite system: the Hilbert space is 
partitioned so that all sites to the left of some bond are in one 
subsystem and all sites to the right are in the other. We find 
one scaling law immediately after the sweep and another in 
the time dependence of the swept state under the final Hamil- 
tonian, in addition to oscillatory behavior. For critical points 
with conformal invariance (z = 1) in one dimension (d = 1), 
the entanglement entropy of a ground state with large finite 
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correlation length diverges as S = f log£ + . . ., where c 
is the central charge of the critical point II 1811 . To obtain the 
entanglement entropy, we have to find the correlation length 
of the state \ip' (g t ,T)) . The process of sweeping generates 
a gap in the system A'(<? t ,r) different from A(g t ) just as 
\if>' Q (gt,T)) is different from IV'u (<?*))■ Polkovnikov has cal- 
culated the scaling form of the "typical gap" A'(g t , T), which 
is A' ~ pWC^+l), Now, combining this with the correlation 
length £' ~ A'<- 2 ) and using z = 1, we obtain 



6(1/ + 1) 



log r + const. 



(3) 



For the transverse Ising model with central charge c = 1/2, 
we find that S = — 1/24 ln(T) + const, which is consistent 
with the overall slope of the numerically calculated entangle- 
ment entropy for the transverse Ising model (</> = 0) just after 
the sweep in FIG.|3ja). A similar expression for the transverse 
Ising model can also be obtained for the entanglement entropy 
of a finite block embedded in an infinite chain lfl9ll . 

If the swept state continues to evolve in the (now con- 
stant) final Hamiltonian, the entanglement entropy oscillates 
around a linearly increasing mean (FIG. Ob)); this linear in- 
crease is as p redicted by Calabrese and Cardy for a "global 
quench" M20I1 . but in the swept case the slope depends on the 
number of excitations created during the sweep. In our numer- 
ical study of the weakly nonintegrable model with <fr = 7r/32, 
entanglement initially grows linearly with time as in the inte- 
grable case but then, after a time related to the interactions be- 
tween "excitations" of the integrable model, rapidly increases 
and makes convergence much more difficult. 

As a consequence of the fact that state \ip' (g t ,r)) includes 
excited states, both S and the expectation values of opera- 
tors that do not commute with the Hamiltonian will oscillate 
as a function of the sweep rate and time. While the sys- 
tem is on the same side of the critical point as tfe, the state 
will evolve adiabatically. Once the critical point is reached 
(at time t c ), the expectation value of any observable that 
mixes the ground state with the band of excited states oscil- 
lates with phase ^jp- w A(g t ). Using the scaling of the gap 
A(fft) ~ \gt - 9c\ zv = T zv {t - t c ) z \ and integrating with 
respect to time, we obtain 
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These findings are consistent with the numerical simula- 
tions of the sweep in the transverse Ising model in FIG. [3] 
The entanglement entropy S shows oscillations in real time 
(FIG. 0b)) which have a constant frequency equal to the gap 
A(g). The oscillations in the rate dependence in FIG. [2a) are 
consistent with Eqn. (O. 

In order to study these oscillations and their damping with- 
out reference to a particular observable, we use a version of 
the Loschmidt echo (see Ref. |2l| and references therein). We 
square the inner product between the wavefunction immedi- 
ately after the sweep and the wavefunction after an additional 




FIG. 4: Exponent a in the overlap \{ip'o(gf,T)\ij}' (gf,T,t))\ = 
exp(al/). The overlap is taken of the wavefunction \ip' (gf, F)} (im- 
mediately after the ramping from g = Qi to g = g j with the rate pa- 
rameter T) and the wavefunction IV'oG?/) T, t)) (after an additional 
evolution for the time t at fixed g — gj). The dots in the upper 
((/> = 0) panel show the analytical results using Eqn. ©. 



wait period of time t: 

l(^o(2/,rM( 5/ ,r,t))| 2 = |<e- iH '% i(9/ ,r)| 2 - e~ a ^ L . 

(5) 

Here we have defined a(t) because this overlap typically 
goes to zero exponentially for translation-invariant states in 
the thermodynamic limit (i.e., when the length of the chain 
L — > oo) unless \tp' (gf,T)) and \ip' (gf,T,t)} are identi- 
cal. This overlap gives a direct probe of the magnitude and 
decay of quantum oscillations in the many-body state and it 
is suitable for computation using iTEBD Il5ll . It has been 
shown to determine the statistics of work done in a quan- 
tum quench lf22Tl . Our main results are that the overlap os- 
cillations have an unusual cusp structure, arising from the 
integrability of the transverse Ising model, and decay even 
in the exact zero-temperature evolution of an infinite system 
(because there is a continuum of excitation energies) alge- 
braically rather than exponentially. 

An analytical approach can be developed from independent 
Landau-Zener tunneling 11231 12411 at each k IU2I1 . In terms of 
the Bogoliubov excitations (see Ref. [12 for details), the wave- 
function \ip'o(gf, T)) for the transverse Ising model (<fi = 0) is 
a product \ifj' (gf,T)) = H k (u k \0) +v k \k,-k)), where |0) 
is the vacuum and the state \k, —k) contains a pair of quasi- 
particles with pseudomomenta (fc, —k). The Landau-Zener 
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formula gives, for the small k's that dominate the tunneling 
for slow rate, 



\v k \=Pk = 1- \ukV =exp - 



2nJ 2 k 2 



r 



(6) 



During the wait period, each fixed-k wavefunction |fc, — k) has 
an energy Af(k) = 2\J A 2 + 4:J 2 gk 2 (this is an approxima- 
tion for low energies), with lattice spacing a = 1. So up to an 
overall phase 

\<(9f,r,t)) = n(^l°) + e " A/(fe) ^l fc '- fc )) ( ? ) 

fc 

Now, using Eqn. (O, the squared overlap in Eqn. (O is rewrit- 
ten as 



n 



dn2 [ A/I*)* , Pk{Pk 1} 
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In the continuum limit, the logarithm of this product becomes 
an integral for a(t). Taking the momentum cutoff to oo, 



1 f c 



dk log 
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A f (k)t 



Pk(Pk - 1) 



(9) 

This is compared to the numerical results in FIG. Ufa). The 
cusplike minima arise from the momentum value k* where 
the tunneling probability P^, = 1/2. These give a singu- 
larity in a(t) at times Af(k*)t = (2n + 1)tt. This sin- 
gular behavior is smeared out in the nonintegrable model at 
4> = it/ '32, which does not have sharp excitations at finite 
energy (FIG. HJb)). This gives a probe of integrable ver- 
sus nonintegrable behavior that could be useful for experi- 
ments [25]. Asymptotics of Eqn. (0 lead to the following 
results: if we have FN <C A 2 for some number of periods 
N > 1, then the overlap has peaks shifted from the points 



At = ttN, with maxima ajv 
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At long times, a approaches a constant independent of the fi- 
nal gap: a(oo) w -0.0564 VF/ J. 

In the integrable case, the amplitude of oscillations in the 
Loschmidt echo falls off as 1/t for large times, which is differ- 
ent from the case of a quasithermal initial state whose excita- 
tion probabilities are given by Boltzmann weights. The latter 
has no cusps and damping as \j\fl. This power-law decay of 
the oscillations, rather than an exponential decay, is an impor- 
tant difference between the evolution found here and the decay 
of some observables after a sudden quench j^]. We note that 
the (pure) state as t — > oo does not resemble a (mixed) ther- 
mal state at any effective temperature of the final Hamiltonian, 
since the excited state occupation probabilities remain deter- 
mined by their tunneling probability Pk, while their energy is 
Af(k); the final state has more occupation of low-fc excited 
states than a thermal state. Sweeping can also induce spatial 
oscillations in correlation functions of local operators 026I1 . 

These results are modified in several ways for non- 
integrable systems. In any system, the Loschmidt echo de- 
pends on both the initial state and the Fourier transform of 



the energy spectrum. The latter is known to change between 
integrable and non-integrable systems: one sharp difference 
in spectral properties is revealed through the smearing out 
of cusps noted above (FIG. |4fb)). More generally, the mea- 
surement of a Loschmidt echo could be used to probe the dif- 
ference in energy level statistics between integrable (Poisson) 
and non-integrable (Wigner-Dyson) systems. 

Local observables are also likely to behave differently in 
integrable and non-integrable systems, e.g., decay of local 
observables to a thermal distribution is expected in the non- 
integrable case. The dynamical crossover to non-integrable 
behavior was observed here as a rapid increase of entangle- 
ment entropy and the destruction of cusps in the Loschmidt 
echo. Accessing the long-time regime and approach to equi- 
librium in the non-integrable case is a major challenge for nu- 
merical methods. 
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